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A FIRM UPPER LIMIT TO THE RADIUS OF THE NEUTRON STAR IN 

SAX J1808.4-3658 
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ABSTRACT 

We show that observations of X-ray pulsing from SAX J1808.4-3658 place a firm 
upper limit of 13.8m^/^ km on the radius of the neutron star, where m is its mass in solar 
units. The limit is independent of distance or assumptions about the magnetospheric 
geometry, and could be significantly tightened by observations of the pulsations in the 
near future. We discuss the implications for the equation of state and the possible 
neutron star mass. 


Subject headings: Subject headings: accretion, accretion discs — binaries: close — 
stars: neutron — equation of state — X-rays: stars 


1. INTRODUCTION 

The discovery (Wijnands and van der Klis, 1998) of an accreting neutron star with a 2.49 ms spin 
period in the soft X-ray transient SAX J1808.4-3658 (orbital period = 2.01 hr, Chakrabarty and 
Morgan, 1998) has important consequences for many aspects of our understanding of low-mass 
X-ray binaries and their evolution. In this Letter we focus on just one of these, namely the rather 
tight limit one can place on the neutron star radius. This in turn has significant implications for 
the equation of state and the possible mass of this star. 


2. THE RADIUS LIMIT 

Our radius limit comes from considering the fact that pulsed X-ray emission is detected at a range 
of X-ray count rates. We make a number of assumptions, which are all commonly adopted by other 
authors; later in the paper we shall consider the consequences of relaxing them. To specify them 
we define the neutron star radius R*, the corotation radius 

Rco = 1.5 X cm (1) 


^Dipartimento di Fisica ed Astronomia, Universita di Palermo, via Archirafi 36, 90123 Palermo, Italy 
^Astronomy Group, University of Leicester, Leicester LEI 7RH, U.K. 



- 2 - 


and the magnetospheric radius 

Rm = 1-9 X cm, (2) 

(e.^. Burderi et al, 1998). Here m is the neutron star mass in solar masses, P 3 is the neutron star 
spin period in milliseconds, c;/!) ~ 1 is the ratio between the magnetospheric radius and the Alfv& 
radius {e.g. Burderi et al., 1998), /U 26 is the neutron star magnetic moment in units of 10^® Gauss 
cm^, Rq is the neutron star radius in units of 10® cm, M 17 is the accretion rate in units of 10^^ g/s. 

We derive a limit using the assumptions Al - A4 defined below. These are that X-ray pulsing 
requires 

Al: 

i?* < Rma.g (3) 


and 

A2: 

.^mag Rcoj (4) 

and in particular that there is no ‘intrinsic’ pulse mechanism (as e.g. in rotationally powered radio 
pulsars). 

We assume also that 

A3: In eqn (^), cj) is independent of M, 

and 

A4: the accretion rate M is strictly proportional to the X-ray count rate F (3-25 keV band) in 
the Proportional Counter Array (PCA) instrument on RXTE. 

Using A4 and A3 we can simplify eqn (|2|) to 

= AF-2/7, (5) 

where A is a constant for a given system. 

X-ray pulsing was detected from SAX J1808.4-3658 at count rates corresponding to Fi = 238 x 
10“^^ erg/cm^/s and F 2 = 21 x 10“^^ erg/cm^/s (Gilfanov et al. 1998). Applying Al, A2 and (|^ 
to these observations gives 

F* < (6) 

and 

AFs"'/" < Fco (7) 

_ 2/7 _ 2/7 

But AF 2 ' = 2.00 X AFr^‘, so by combining the two we get the extremely simple result 

F* < 0.5iico = 13.8 X km ( 8 ) 
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Usually our lack of knowledge of the uncertain quantity (p and to some extent M frustrates attempts 
to get precise quantitative information from assumptions like Al, A2. However here all of these 
uncertainties are removed because of the scaling (|5|): given the assumptions Al - A4, the upper 
limit (HI) depends only on the neutron star mass. 


3. IMPLICATIONS 


Neutron-star equations of state generally yield a theoretical mass-radius relation of the form 

(9) 

When combined with such a relation the limit (^) defines a minimum mass for the neutron star in 
SAX J1808.4-3658 , while the characteristic lengthscale L measured from the Type I X-ray bursts 
(in’t Zand et al., 1998) defines maximum mass via the condition L < (i.e. the burst site cannot 

exceed the entire neutron-star surface). This last condition is derived assuming that the peak 
luminosity during the X-ray bursts is close to the Eddington luminosity: 


L = 31.3 X (kT/lkeV) km < R^ 


( 10 ) 


In applying equation (^) two points have to be considered. Firstly, the blackbody temperature 
(fcT = 2.23 keV) was determined by averaging the values obtained from the first 7 sec of the burst 
reported in’t Zand et al. (1998) together with the first 7 sec of the second burst. During the initial 
phase of the bursts photospheric expansion is likely, so the photospheric radius L could exceed 


the radius of the neutron star. Secondly (H) is derived under the assumption of pure blackbody 
emission, while more accurate determinations of the spectrum of the emitting photosphere show 
that several corrections apply to the spectral shape (see e.g. Titarchuk, 1994), resulting into a 
underestimate of the radius of the photosphere. Figure 1 shows a mass-radius plane with these 
constraints marked for various equations of state. For any reasonable equation of state the neutron 
star radius must be smaller than 15 km. As can be seen, very stiff equations of state (EOS) (such 
as the Mean Field EOS labeled as L see e.g. Cook, Shapiro and Teukolsky, 1994) are marginally 
allowed if the m takes the favoured value of 1.4, and softer equations of state (such as the Reid soft 
core EOS, labeled A or the FPS EOS see e.g. Cook, Shapiro and Teukolsky, 1994) can currently 
accomodate even smaller masses. However the upper radius limit (H|) is at present only derived for 
X-ray count rates down to T 2 = 21 x 10“^^ erg/cm^/s. The system is detected at much lower rates, 
and it will be possible to check for pulsation and thus reduce the allowed area in the R* — m plane 
quite drastically in the very near future. The Figure shows the resulting constraints if pulsing is still 
detectable at count rates 9 x 10“^^ erg/cm^/s and 3 x 10“^^ erg/cm^/s. The detection of pulsing 
at a count rate 3 x 10“^^ erg/cm^/s would force M to be significantly larger than 1.4 for even the 
softer equations of state. This would be quite reasonable, given that the formation mass for the 
secondary star in neutron-star LMXBs is probably IMq (cf King &: Kolb, 1997, Kalogera, Kolb 
&: King, 1998), and that significant mass must have been transferred to the neutron star during its 
evolution. 
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4. DISCUSSION 


Here we consider the effect of dropping each of the assumptions A1 - A4. 

A1 is normally regarded as the minimum condition for the magnetic field to be dynamically im¬ 
portant. Pulsing while this condition was violated would suggest that the neutron star had some 
intrinsic pulse mechanism. This is very unlikely to be the usual radio-pulsar process, as this is 
known to be effectively extinguished if there is significant matter inside the light cylinder, as would 
be true in this case. Moreover the idea of such an intrinsic component taking over from the normal 
accretion power-pulsing at low count rates is belied by the almost constancy of the pulse fraction 
(Gilfanov et al, 1998) 

Dropping A2 would imply that pulsing was possible in the ‘propeller’ regime, when matter is 
thought to be centrifugally expelled rather than accreted. Again it is hard to reconcile this idea 
with the constancy of the observed X-ray pulse fraction. 


Assumption A3 appears to be perhaps surprisingly accurate, considering that it is only a dimen¬ 
sional estimate that is given by comparing magnetic and material stresses. Campbell (1997) derives 
a general prescription for the disruption radius i?mag,disc of a thin disc of central density p by a 
dipole aligned with the spin axis 

^mag.disc OC (11) 


The weak p-dependence in this relation means that if disruption occurs in a region of the disc 
where gas pressure is dominant and thus p oc (cf Frank et ah, 1992) one gets i?mag,disc oc 

i.e. quite close to the dimensional estimate (|^ with cj) ~ constant. If instead radiation 
pressure is important near the disruption radius one has p M which would give the 

stronger dependence 


R 


mag,disc 




( 12 ) 


However the rather low luminosity (^ 10^® erg s“^) of SAX J1808.4-3658 implies accretion rates 
M ^ 10^® g s“^, so that radiation pressure only becomes comparable with gas pressure for disc radii 
R ^ 20 km (Frank et al., 1992), as compared with the likely corotation radius Rco — 30 km. The 
disruption therefore probably does not occur in a region dominated by radiation pressure. However 
for comparison we plot the curve for F 2 = 21 x 10“^^ erg/cm^/s with the dependence (12) in Figure 
1. As is possible to observe, this curve is well below the lower limit imposed by (0), supporting 
our hypothesys that the disc is truncated in a zone where radiation pressure is negligible. 


The X-ray spectrum of SAX J1808.4-3658 is remarkably constant (Gilfanov et al., 1998), suggesting 
that A4 is unlikely to be upset by a changing bolometric correction. A more serious possibility 
is that some of the accretion energy might be stored in other forms, rather than being released 
promptly as radiation. For example Priedhorsky (1986) suggested that the neutron star spin might 
provide such a reservoir in the QPO sources. To weaken the limit (^) significantly would require 
essentially that the drop in X-ray count rate masked an almost constant accretion rate. Again the 
constancy of both the X-ray spectrum and the pulse fraction tend to go against this idea. 
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We conclude that the limit ^ is robust. The work of this paper also shows that the magnetospheric 
radius in SAX J1808.4-3658 lies between 1 and at most a few times ii*- While it is possible that 
this is merely a coincidence, it may instead suggest that whatever process induces field decay tends 
to stop once such conditions are achieved. This in turn supports the idea that the decay may be 
induced by accretion. 

After submitting our paper for publication we became aware of an independent paper by Dimitrios 
Psaltis and Deepto Chakrabarty which reaches similar conclusions. 
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FIGURE CAPTIONS 


Constraints on the radius of the neutron star in SAX J1808.4-3658 . This is bounded above by the 
solid curve labelled 7.0 x 10^® erg s“^, corresponding to the lowest X-ray luminosity at which pulsing 
is currently detected (and derived adopting the assumption A3, see text), and bounded below by the 
size derived from the observation of Type I X-ray bursts (lower solid line). Detections of pulsing at 
3.0 X 10^® erg s“^ or 1.0 x 10^^ erg s“^ would reduce the upper limit to that given by the appropriate 
dashed curve. For comparison, the upper limit derived adopting the disruption radius dependence 
([T^ and the luminosity 7.0 x 10^® erg s“^ is shown by the curve labelled Radiation pressure limit. 
The curves ‘L’, ‘FPS’ and ‘A’ give the mass-radius relations predicted by the equations of state 
discussed in the text. For each such equation, the upper radius limit defines a maximum radius 
and minimum mass for the neutron star. Thus equation of state L currently requires M ^ 1.34Mq, 
R ^ lb km. The latter is probably the maximum reasonable radius for the neutron star in SAX 
J1808.4-3658 . 
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